Large iV reduction in the continuum three dimensional Yang-Mills theory 
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Numerical and theoretical evidence leads us to propose the following: Three dimensional Euclidean 
Yang-Mills theory in the planar limit undergoes a phase transition on a torus of side I = l c . For 
I > lc the planar limit is /-independent, as expected of a non-interacting string theory. We expect 
the situation in four dimensions to be similar. 



Introduction. Yang-Mills theory in three dimensions 
is similar to Yang-Mills theory in four dimensions in ex- 
hibiting a positive mass gap, linear confinement, a fi- 
nite temperature deconfinement transition and a sensible 
large N limit 0]. 

Nonabelian gauge theories in 3 and 4 dimensions inter- 
act strongly at large distances and weakly at short dis- 
tances. Doing calculations that bridge these two regimes 
remains a major challenge, of central importance to par- 
ticle physics. It has been a long held hope that the task 
would simplify at infinite number of colors, N. Here, at 
strong coupling, a fitting hypothesis is that the theory 
describes free strings, while at weak coupling, the theory 
certainly describes weakly interacting particles. The re- 
sult of this paper indicates that strong and weak coupling 
are separated by a phase transition at infinite TV. Specif- 
ically, we provide numerical evidence that SU (N) gauge 
theory on a three torus of side I undergoes a transition at 
a critical length, I = l c . For any finite value of N there 
cannot be any phase transitions in this system. The ex- 
istence of such a transition at infinite N is surprising, 
raises questions about the usually assumed smooth de- 
pendence of observables on momenta and might indicate 
deeper connections between gauge theory, string theory 
and random matrix theories. 

Over twenty years ago, in the context of SU(N) lattice 
gauge theory, Eguchi and Kawai made the observa- 
tion that at infinite number of colors space-time can be 
replaced by a single point. This dramatic reduction in 
the number of degrees of freedom should make it easier 
to deal numerically with planar QCD than with ordinary, 
three color QCD. For a practical procedure it is essential 
that some version of large TV reduction also hold in the 
continuum, not just on the lattice. Previous attempts 
to define a continuous reduced model had problems with 
topology and fermions [3j. 

We focus on the continuum limit of pure lattice YM 
defined on a torus and try to establish that expectation 
values of traces of Wilson loop operators do not depend 
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on the size of the torus. Wilson loops of arbitrary size can 
be folded up into the torus and correctly reproduced p| . 
The lattice is essential because it provides a regulariza- 
tion with well defined loop equations Q . Loop equations 
provide a convenient tool to establish reduction. 

We restrict ourselves to three dimensional theories for 
numerical reasons. We find that continuum large N re- 
duction holds so long as the torus is large enough. The 
critical side length of a symmetrical torus is denoted by 
l c and is defined in terms of a microscopic fundamental 
physical scale of the theory. Solving the theory for some 
I > l c would produce complete and exact information at 
leading order in N for any I. The system as a whole 
undergoes a phase transition at / = l c . The number of 
sites in a numerical simulation in a given direction, L, de- 
termines the maximal value the ultraviolet cutoff A can 
take. It is A = i ± . For a Wilson action the lowest L 

'c 

that has some semblance to continuum is L = 3. Thus, 
at the expense of larger N one can get numerically close 
to continuum using very small lattices. The values of N 
needed are of order 20 to 50 and this trade-off is worth 
taking. 

If a similar result holds in four dimensions, a shortcut 
to the planar limit becomes a realistic option. Our expe- 
rience makes us hopeful and our tools should allow us to 
tackle four dimensions in the future. 

A lattice argument. There is a global Z d (N) {U d {\) 
in the N — > oo limit) symmetry on the torus that leaves 
contractible Wilson loop operators invariant but multi- 
plies Polyakov loops winding around a direction /i by a 
phase e~ k >* . The preservation of this symmetry is cru- 
cial for large N reduction 0. Eguchi and Kawai have 
shown that the lattice loop equations in the N = oo 
limit on a single site lattice are the same as on an infinite 
four dimensional lattice as long as the U 4 (l) symmetry is 
unbroken. The continuum limit in the single site lattice 
model has to be taken by sending the lattice coupling 
b = t° infinity, but in d > 2 a phase transition 
occurs, blocking the way. At the transition the Z d (N) 
symmetry breaks spontaneously, ruining the equivalence 
of loop equations. It is possible to fix the single site lat- 
tice model by quenching or twisting y\ the system. 
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We take a different approach here. The proof of Eguchi 
and Kawai goes through for a lattice torus of size h xfa— 
with arbitrary and in any dimension. The loop equa- 
tions, together with boundary conditions for small loops, 
establish equality of expectation values of traces of op- 
erators associated with arbitrary finite closed loops in 
the infinite volume theory and their folded, contractible, 
counterparts on the torus. Suppose we reduced the model 
to only an L d lattice with L > 1: Again we expect the 
global symmetry to break if b > b c (L) and reduction will 
hold for b < b c (L). b c (L) will increase with L and if b c (L) 
depends asymptotically on L as dictated by microscopic 
scaling for d = 3, 4 then continuum large N reduction 
will hold if we take the limit by keeping b < b c (L) and 
taking b — > oo. 

In the approach pursued here, we have to deal with 
one lattice artifact. There will be a cross-over in the 
lattice internal energy for the Wilson gauge action at 
some small b for a finite torus size and a finite N. 
The cross-over becomes a "bulk" transition at infinite 
N, occurring at bf(N = oo, L) for any finite lattice 
of size L d in lattice units. Lattice large N reduction 
would imply that bf(N — oo,L) does not depend on L, 
bc(N — oo, L) = b~ . This is consistent with numer- 
ical simulations. The loop equation, together with con- 
straints which come from the parallel transporters being 
unitary matrices, produce the "bulk" transition without 
loosing their validity or changing their form. The lat- 
tice transition occurs when the unitary matrix associated 
with the one by one loop opens a gap at eigenvalue -1 in 
its spectrum in the large N limit. As b increases fur- 
ther the gap widens. In the continuum this means that 
parallel transport round a tiny loop will not differ much 
(in norm) from the identity. Similar transitions occur at 
bc(N,L — oo) for large enough N. The common limit- 
ing value at bf(N — oo, L = oo) = bf is rapidly ap- 
proached. This family of transitions are lattice artifacts 
not associated with any symmetry breaking. Examples 
are the Gross- Wittenlfj transition in two dimensions and 
Creutz's transitions [9j for N > 4 in four dimensions. 

Even though lattice reduction is valid on either side of 
b~(N — oo, L) as long as one is below b c (L), we have 
to be above b~(N — oo,L) to realize continuum reduc- 
tion. For L = 1 (the Eguchi-Kawai model) and d > 2, 
the infinite N "bulk" and Z d (N) breaking transitions ac- 
cidentally fuse at a b c ^ b^°°. Similar "accidents" can 
happen for L=2,3.., but a window opens for large enough 
L between b^ and b c (L). In three dimensions, an L = 3 
lattice already has a window. 

6f°° = 0.5 .8) and b c (L) = oo in d = 2. The 
U 2 (l) symmetries are not broken and continuum reduc- 
tion works on tori of any size in two dimensions. The 
"bulk" transition occurs close to bf = 0.4 in d = 3. Or- 
dinary scaling in d = 3 would require yTl) ^° a PP roacn 
a finite nonzero limit as L — > oo. Monte Carlo simula- 
tions were performed using a combination of heat-bath 
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FIG. 1: Eigenvalue density distribution of a 4 x 4 Wilson 
loop on 4 3 (folded) and 6 3 (unfolded) at 6 = 0.66 and N = 23. 

updates by SU(2) subgroups and of full SU(N) over- 
relaxation steps. Ultraviolet fluctuations in loop observ- 
ables were suppressed by APE blocking [ljj. We mon- 
itored the eigenvalue distribution of the Polyakov loops 
in the three directions and found that 0.6 < 6 C (3) < 0.7, 
0.8 < 6 C (4) < 0.9, 1.0 < 6 C (5) < 1.2 and 1.2 < 6 C (6) < 
1.35. When combined, these results indicate that the 
scaled critical coupling j-pTj is in the region [4.2,5]. We 
compared folded and unfolded versions of the same loop 
on tori of different sizes and found the spectral densities 
associated with them to match as long as > 5. An 
example of such a comparision is shown in Fig. ^ We 
also checked scaling by comparing Wilson loops of same 
physical size at different lattice spacings. An example of 
scaling is shown in Fig. [21 

Continuum perturbation theory. If we had a 
scalar field theory where the field is a hermitian N x N 
matrix we know that independence on the torus size is 
impossible. This dependence does not go away in the 
planar limit. On the level of Feynman diagrams (taken 
in coordinate space) it is easy to calculate the depen- 
dence on torus size for large /, in particular if the theory 
is massive |l l| . The leading correction is exponentially 
suppressed in I and comes from one virtual particle going 
round a non-contractible circle on the torus. There is a 
stable particle like that and it is in the adjoint representa- 
tion of SU (N). In the gauge case, if there is confinement, 
we could use only singlets under SU (N) and, at infinite 
N, these singlets make sub-leading contributions to the 
free energy at leading, 0(N 2 ), order. We conclude that 
for a confining gauge theory a planar diagram with a rib- 
bon (double- line) representation of propagators makes no 
contribution if one tears one of the propagators out of the 
surface and winds it round the torus. 

Another way to see how reduction works in perturba- 
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FIG. 2: Eigenvalue density distribution of L x L Wilson loop 
on L 3 for L/b = 5 and L = 4,6. JV is set to 23. 



tion theory is to understand what happens to momen- 
tum space @. Having a torus means that momenta are 
quantized in units of =S and there is no way around 
this for a massive scalar matrix field. In the gauge case 
the Feynman expansion starts from a constant gauge 
field background. The gauge invariant content of this 
moduli space consists of d sets of angles (9* which ef- 
fectively fill the intervals between the quantized mo- 
menta making momentum space continuous and I inde- 
pendent. The filling has to be uniform and this is true 
at infinite N if the global Z d (N) symmetry is unbro- 
ken. The background - in a translation invariant gauge 
- is given by = diag{9 1 9 2 , ....O 1 ^) but only the set 
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^ } labels truly distinct vacua. This is 
9 l \ < j, exactly as needed to fill in the gaps be- 
tween the At infinite N the vacuum is character- 
ized by d eigenvalue distributions on the complex unit 
circle. 

Hints from string theory. In view of developments 
during the last few years [12j h seems more likely now 
than ever before that indeed large N SU (N) pure gauge 
theories are equivalent to some string theory at zero 
string coupling. This means that the logarithm of the 
partition function defined on a finite torus and divided by 
the volume of the torus, is, in the planar limit, given by a 
sum of extended, spherical, two dimensional excitations 
embedded in the same torus. But, there is no way for the 
spherical surface to become non-contractible on the torus 
and thus it cannot detect that target space is a torus 0] . 
Hence, one can have no dependence on I. It is well known 
that simple string models on toroidal backgrounds can- 
not distinguish very large radii from very small ones; l c , 
as a minimal radius, realizes a similar phenomenon in the 
unknown non-interacting string theory describing planar 
three dimensional pure YM. 



It used to be revolutionary to think that statistical 
field theories on finite volumes can have phase transi- 
tions. This is no longer true. To the early toy model 
examples |3 we can add now cases of true, full fledged 
field theories with real relativistic degrees of freedom, also 
developing phases transitions in the planar limit |12| . 

Large N phase transitions. Large N transitions 
may emerge as quite ubiquitous in continuum gauge the- 
ories. There are transitions, like the one presented in this 
paper, that affect the system as a whole, but there are 
also other transitions that affect only a class of observ- 
ables [3ESE1E3- 

The basic observables used in our study have been the 
distribution of eigenvalues of Wilson and Polyakov loops. 
For large N these observables are unconventional because 
they involve traces of all powers of the basic unitary ma- 
trix, not only a few low powers. Thus, issues of renor- 
malization require more work [J^. If these issues can 
be resolved, we might be able to exploit the fact that 
eigenvalues of large matrices have many universal prop- 
erties [2(|: The dynamics of the gauge theory could be 
encoded in the transformations one needs to carry out 
in order to bring these eigenvalue distributions to uni- 
versal forms. While there are difficulties in continuum 
perturbation theory, the situation on the lattice is very 
clear: We numerically look for features that scale as the 
universal features of the field theory would have it. 

The simplest strong-weak transition would be associ- 
ated with Wilson loops: Small loops will have paral- 
lel transporters with a spectral gap and big loops will 
have almost uniform distributions to account for confine- 
ment in all finite irreducible representations. At finite N 
there are no gaps in the spectra but, in the range of the 
would be gaps, the eigenvalue density is exponentially 
suppressed as N increases. 

Beyond the transition. For I just a bit smaller than 
l c exactly one of the Z(N) factors in the Z 3 (N) breaks 
spontaneously. Thus, the forty eight clement cubic sym- 
metry group of our equal sided torus breaks down to an 
eight dimensional group acting in the plane perpendicu- 
lar to the direction in which the Polyakov loop spectra 
took on non-uniform structure. 

In order to prepare ourselves for what to look for when 
the torus is further squeezed we studied the l d EK model, 
now interpreted as a simple effective model for the dy- 
namics of the vacuum manifold of the full system. Simu- 
lations we have carried out in three and four dimensions 
showed that at infinite N these models undergo a stair- 
case of transitions, breaking one additional Z(N) factor 
at a time. The possible continuum meaning of the various 
intermediate phases will have to wait for more work. 

In super-symmetric YM gauge theories, compactificd 
super-symmetrically on tori, the perturbative mechanism 
driving the spontaneous breaking of the Z d {N) symmetry 
can be eliminated. Beyond perturbation theory we do 
not know the answer, and other global symmetries come 
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into play. It is conceivable that in some cases l c = 0, 
indicative of a pure matrix model representation of the 
planar limit of a continuum gauge theory. Although the 
physical size is zero, regularization issues might require 
one to take L — > oo in a way correlated with b — ► oo, 
and the zero size model may not admit a definition as 
the large N limit of an ordinary matrix integral. 

Future lattice work. Building on earlier two dimen- 
sional work we know how to calculate meson propagators 
in the planar limit. Meson momenta of values below the 
ultraviolet cutoff can be introduced by multiplying the 
original link matrices U^{x) by phase factors e lJV . The 
p M allow to tune the momenta carried by the mesons to 
desired values. One has no finite volume effects to worry 
about: to get to the continuum limit one just increases 
b, making sure that I stays larger than l c . Values of N in 
the range of few tens seem to be adequate. The lattice 
Dirac matrices are much smaller (and much denser) than 
in usual simulations. We would be able to address the 
smoothness of the two point meson correlation function, 
at infinite N, as a function of q 2 , where q is Euclidean 
momentum. Could there be a non-analyticity at some 
q 2 l After all, if the crossover between physical strong 
and weak gauge forces happens in a range of scales that 
shrinks to zero at infinite N, phase transitions may oc- 
cur in every observable, not only special ones, like Wil- 
son loops. In four dimensions this could, finally, bring 
about a peaceful coexistence between large N and in- 
stantons |21| . 

In parallel simulations of the pure gauge case, large N 
work will require a floating point effort per node that 
grows as N 3 while communication demands will only 
grow as N 2 . So, PC farms with off the shelf commu- 
nications would be well suited. 

Conclusions. Our final conjecture about three di- 
mensions is stated in the abstract. We call it a "con- 
jecture" because our numerical tests have been relatively 
modest and because the consequences of the conjecture 
could be far reaching: many applications of 't Hooft's 
large N limit [22] assume analyticity in momenta and 
this assumption is now challenged. Our evidence is a 
combination of numerical work and more theoretical ob- 
servations. On the numerical side we see the Z 3 (N) sym- 
metry breaking point on the lattice change with lattice 
size in a way consistent with continuum scaling. Theo- 
retically, lattice large N reduction based on large N loop 
equations is a strong coupling argument while averaging 
over the moduli space of constant abelian connections at 
weak coupling resolves an apparent contradiction with 
conventional wisdom about finite size effects. 
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